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Abstract 

We consider free electrons in rectangular quantum dots, with either hard wall boundary condi- 
tions or anharmonic confinement. In both cases, due to finite size effects, a homogeneous electric 
field applied along one of the rectangular axis is shown to induce abrupt changes in the electron 
density, parallel and perpendicularly to the field direction: the electron density jumps from one 
configuration to another. Making use of this property, we propose a purely electrical mechanism to 
control the magnitude of the effective exchange coupling between two quantum dots. This system 
has been proposed recently as a quantum gate for quantum computation. 

PACS numbers: 73.22.-f, 85.35.Be, 03.67.Lx 
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I. INTRODUCTION 



Due to remarkable advances in microfabrication techniques, it is possible nowadays to 
study the spectroscopic properties of single systems at mesoscopic or even nanoscopic scale. 
These can be semiconductor quantum dots,^^ ultrasmall metallic particles^ or, organic 
molecules,— for instance. This opens many new perspectives for fundamental physics but 
also for the development of new technologies. Indeed, on the one hand, it is possible to study 
properties of very small systems where quantum phenomena are prominent; on the other 
hand, one may think of applying these new properties to go one step further on the way 
to extreme miniaturization towards electronics at the nanometer scale. In this context, to 
understand the behavior of small systems under the influence of external excitations is a cen- 
tral issue, with the aim to control the parameters of future nanodevices by external sources. 
Here we focus our attention to DC homogeneous electric fields applied to semiconductor 
quantum dots. Since in all the above mentioned systems electrons are confined in a very 
small region of space giving rise to a discrete electronic spectrum, we expect they should 
share common features, despite important differences. Therefore, we believe some of the 
conclusions described below are relevant for metallic particles and large organic molecules, 
as well. 

The problem of the influence of a DC electric field on electronic systems has attracted a lot 
of interest since many years. In any case, we don't want to review this tremendous amount 
of articles but, instead, we content ourselves to cite those of direct use for us. Most often 
in the literature, this problem is treated in one- dimension, either for free electrons^i^ or for 
electrons in a periodic potential,^^ with the main focus on the Wannier-Stark ladder problem 
in the latter case. Here we investigate systems with many channels. More precisely, we don't 
consider Zener tunneling between bands, but we are rather interested in effects caused by 
the interplay of the applied electric field that distorts the electronic bands and the Pauli 
principle which acts kind of quantum pressure. 

We consider rectangular semiconductor quantum dots with characteristic sizes such that 
their mean level spacing ranges in the meV region. The homogeneous electric field is applied 
parallel to the long axis. The electrons within the box are supposed to be non-interacting; in 
particular, the electric field is assumed to be unscreened. This is a good approximation in the 
case of strictly one dimensional quantum dots, as it was shown recently for molecular wires.— 
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It should also be valid for quasi-one dimensional dots where only few channels are occupied 
which is precisely the case of interest in the present work. A first obvious effect of the 
field is to polarize the electronic cloud in the longitudinal direction i.e. parallel to the field 
direction. Our purpose is to investigate another effect caused by the finite size of the system. 
We consider two types of confinement: infinite potential wall and harmonic confinement plus 
a small anharmonic perturbation. In both cases, the electric field is shown to induce at a 
critical value, Ec, abrupt changes in the electron density in the two directions, parallel and 
perpendicular to the field. Despite the fact that no interband coupling exists, the electric 
field modifies the distribution of the electrons among the different bands. The net effect 
of this redistribution is to decrease the longitudinal polarization and to push the electrons 
toward the surfaces of the well perpendicularly to the field. This effect is discussed here on 
the basis of simple models, with the purpose of pointing out the basic mechanism at its origin 
and its possible relevance in the context of nanophysics and quantum computation. Indeed, 
as a first application, we consider two lateral weakly coupled quantum dots: the electric 
field is shown to modify abruptly the gap between the lowest singlet and triplet states. This 
last effect could have important application in the context of quantum computation were 
similar devices where proposed as possible realization of quantum gates. ^ 

The paper is organized as follows. In section II we investigate the case with hard wall 
boundary conditions, and in section III the case with harmonic confinement plus a perturba- 
tive anharmonic term. Finally, in section IV, we present our results for two weakly coupled 
quantum dots before closing with some conclusions. 

II. RECTANGULAR POTENTIAL WELL WITH INFINITE WALLS IN AN AP- 
PLIED ELECTRIC FIELD 

We start by considering a system of free electrons confined in a rectangular potential well 
with infinite walls. The corresponding one-electron time-independent Schrodinger equation 
in presence of an applied electric field is 



where m is the effective mass of the quasi-electrons, q the electron charge (g < 0) and 
E the applied electric field. The wave functions fulfill the hard-wall boundary conditions: 




(1) 
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\I/(0, y) = "^{Lx, y) = and 0) = Ly) = 0. We assume a rectangular box such that 
> Ly (cf. Fig. 1). 

Without electric field {E = 0), we recover one of the basic problems of quantum mechan- 
ics: the well known particle in the box problem in two dimensions. The eigenenergies are 
given hy^ 

W. = |^[(^r + (^r] (2) 

X y 

where {nx,ny} G N* x N* (N* is the set of integers without 0). 

With electric field {E ^ 0), since there is still no coupling between the x and y coordinates, 
the wave function may again be written as a product: "^{x, y) = ip{x)C,{y)- The equation for 
C,{y) is the one of a free-particle in a one dimensional box.^^ The equation for (f{x) reads 

;^{x) + qEx(p{x) = e(p{x). (3) 



2m dx"^ 

Solutions of this problem are well known for both infinite^ and finite size systems.-^ We 
closely follow these two references. 

We set e = and F = —^qE and write down the solution of Eq. Q as a linear 
combination of Airy functions,— Ai and Bi, 

V{x) = «Ai( - F'/'[x + ^]) + (3Bi{ - FV3[a; + |,]) (4) 

where a and (3 are two real constants to be determined. The wave functions should vanish 
at the frontiers of the potential well, V'(O) = ^{Lx) = 0. These conditions help on the one 
hand to fix the values of a and P, and on the other hand to write down the secular equation 
that gives the set of discrete eigenvalues, e, 

Ai( - gF-2/3)Bi( - _ _ _ L^^i/s _ eF-^/3)Bi{ - eF-^'^) = 0. (5) 

We can equivalently express the Airy functions in terms of Bessel functions of fractional 
order,— J1/3 and J_i/3, 

Ai(-^) = f (Jl/3(C) + J-1/3(C)) 
Bi(-^) = ^(J_1/3(C)-Jl/3(C)) 



with C = f-z^''^- The eigenvalue equation can then be rewritten as 

(7) 



The equations (jS)) and ((Tj) can be solved numerically, but additional approximations may 
be applied for our purpose. Indeed, especially at weak electric field, most of the eigenvalues 
are such that e'^'^ ^ FLn,^. In this case, we can make use of an asymptotic expansion of the 



Bessel functionAl^ (when z — >■ +cxd, ~ ^J2/txz cos(z — yU7r/2 — 7r/4)) that greatly simplifies 
Eq. dZj) to yield 

[L,F^'^ + - ^'^ = (8) 

with n G N*. The latter equation is easier to solve than Eq. ((Tj) and is very reliable for 
large n. In fact, for all the cases we have considered, the results obtained from Eq. (jH)) are 
in excellent agreement with the exact ones, even for small n. To give an example, we get 
for FL;j. = 10~^ a relative error of 10~^ for the first eigenenergy, n = 1. The error rapidly 
decreases with increasing values of n. Starting from Eq. we may further assume that 
e ^ FLx- We can then expand the first term of Eq. (jHJ in powers of At first order 
one gets, as expected, the free electron spectrum e ~ (f^)^- Details about the opposite 
limit where e <^ FL^. can be found in Ref. 0. In any case, the spectrum of the system with 
electric field is given by 

Examples are shown in Fig. 2, for L^. = lOLy. The two lowest bands, Uy = 1 and Uy = 2, 
are shown with and without electric field. One can notice that the distortion of the bands 
induced by the electric field, depends on the energy or quantum number n^: it is rather 
strong at low energies (small n^) where a change of curvature can be seen in our example, 
and turns into a uniform shift at higher energies (large n^.). 

At zero temperature, the ground state of the system with Ne electrons is determined 
by filling up successively the one-electron eigenstates according to the Pauli principle. In 
analogy with classifications done in atomic physics, it may be characterized by a set of 
occupation numbers, Any{E), that give the number of electrons in band ny at a particular 
electric field. The ground state may then be characterized by the product Yluy An 
meaning that the band Uy is populated by A^y electrons; it gives a representation of the 
ground-state electronic configuration. Since the magnitude of the perturbation is not the 
same for every state f{x), one may expect the electric field to be able to induce changes in 
the set of occupation numbers. Such an example is seen in Fig. 2, where we have considered 
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34 electrons, 17 with spin up and 17 with spin down. Without electric field, only states in 
the lowest band are doubly occupied (this is represented by the open diamonds in Fig 2) 
and in such a way that the lowest unoccupied level is in the second band. With an electric 
field, the perturbation of the lowest level of the second band is larger than the perturbation 
of the highest occupied level of the first band. Consequently, an abrupt change in the set 
of occupation number occurs once the electric field becomes larger than a critical field, Ec, 
for which the lowest occupied and the highest unoccupied levels become degenerate. The 
value of this critical field dramatically depends on the geometry of the box: we obtain for 
instance, Ec ~ ^■^2^^\ ^'-'^ ^ lOLy and Ec ~ ^■^^2^^\ = 9-8Ly. For the case 

shown in Fig. 2, the ground state is characterized by the product l'^'^ for E < Ec, and l'^^2^, 
for E > Ec. 

Since a particular band, characterized by Uy, is associated with a particular function 
^ny{y) — -^/^^^^(t^^)' ^ change in the occupation numbers, corresponds a change in 
the electron density in the direction perpendicular to the applied electric field. This is well 
seen from the transverse electron density, p^ly, E), which is the density with the x variable 
integrated out 

p^{y,E) ^ dxp{x,y,E) = tt^ E (^) T^^ (10) 

Uy 

where the upperscript T is for Transverse. This quantity is expressed by the sum of the 
occupation numbers, An^, weighted by the square of the corresponding transverse wavef unc- 
tion, $^ny. An example is shown in Fig. 3, where we have plotted the change of transverse 
density defined as the difference of E) with and without electric field 

5p^{y, E) = p^{y, E) - p^(y, 0) = ^ K(^) - A„,(0)] sin^ ^-^y. (11) 

^ ^ Uy ^ 

Since Any{E) — A„^(0) for E < E^ this is a discontinuous function of the electric field, non- 
vanishing for E > Ec only. In our example {L^ = WLy, Ne = 34), at E = Ec = two 
electrons jump from the band n^^ = 1 to the band riy = 2 inducing the changes seen in Fig. 3: 
electrons leave the center of the dot to reach the two edges parallel to the applied electric 
field. Of course, this jump produces also a decrease of the longitudinal electric-polarization 
i.e. in the x-direction, an example is given in the next section. 
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III. ANHARMONIC OSCILLATORS WITH AN APPLIED ELECTRIC FIELD 



As a next example, we consider electrons within a two dimensional box as in the previous 
section, but with a more realistic harmonic confinement. We assume the potential to be 
asymmetric (rectangular box) and we add a weak anharmonic term in the direction of the 
applied electric field. The corresponding one-particle time-independent Schrodinger equation 
is 

11 

-— V^^'(x, y) + [-mulx^ + -mujly'^ + Xx'^ + qEx]^!{x, y) = e^(x, y) (12) 

where we assume ujx < ojy and A > 0. In practice, due to irregularities in the confinement 
potential anharmonicity always exists. Another alternative to justify this term in the model 
(fT^ . would be to consider it to have a more fundamental origin, resulting from a renormal- 
ization procedure of a more realistic model including Coulomb interaction, for instance. In 
any case, the anharmonic term introduces a field-dependent perturbation, function of the 
unperturbed energies, which does not exist in the pure harmonic case and which is crucial 
for our purpose, as will be explained in the following two subsections. 



A. Harmonic oscillators in an electric field 



First, we neglect the anharmonic term and consider the case without electric field, A = 
= 0. In this case one recovers the very well known two dimensional harmonic oscillator 
problem. The wave-functions are written as a product, ^nx,ny 

{x,y) = ipnSx)ipny{y)-, with 

^nM = i^)'^'^^e-^'^-'/'H^^iM (13) 
vr v2""n„! 

where u = x/y, = ■\f^^^ the inverse of the effective Bohr radius, and is the Hermite's 
polynomial of order n^. The corresponding eigenenergies are 

e„„ny = {nx + ^)^c^x + {ny + "^^^y (14) 

where Ux and ny are two integers.-^ 

With an applied electric field but without anharmonicity, A = and 7^ 0, the overall 
spectrum of the harmonic oscillator is shifted by — which corresponds to a shift of the 
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center of mass by — in the x-direction. The eierenf unctions become^ 

The harmonic confinement and the shift induced by an homogeneous electric field are 
schematically shown in Fig. 4. Since the perturbation induced by the field is the same 
for all eigenenergies, there is no change in the set of occupation numbers, Any , and thereby, 
in the electron density. This property is very peculiar to the harmonic potential confinement. 
The situation is very different with anharmonicity as it is shown in the next subsection. 



B. Effects of anharmonicity 

We assume A very small compared to the other energy scales, A//5^ <C fioju {u = x/y). It 
has been known that the perturbative expansion in A is not convergent but asymptotic— 
Yet, the perturbative results may approximate the exact one if they are restricted to low 
orders. Here we assume A sufficiently small to consider the anharmonic term at first order in 
perturbation theory only. Although the range of validity of this treatment is very restricted, 
our purpose is mainly qualitative and aims at pointing out basic mechanisms rather than 
reaching a quantitative description of the phenomenon that would require more sophisticated 
approaches.— We believe that the effects describe here occur at any A. We calculate the 
correction 

/+00 r+oo rp 

dx^l{x,E)x' = \ / dx^l{x,0){x-^^Y. (16) 

Exploiting the well known recurrence relation for Hermite polynomials, xHn{x) = 
nHn-i{x) + ^Hn+i{x), and noting that the integrals with odd power of x are null, we get 

, 3 A/i^ , 2 In In 

An. OS = 7:^^{nl + n, + -)+ 6A^Vt^- + o 



plus an overall shift of the spectrum A ^4^3 that plays no role here; thus we neglect it in 
the following. The last term of (fTTj) is driven by anharmonicity and introduces an explicit 
dependence on the applied electric field, and on the quantum number, n^: an applied 
electric field is able to change the electronic configurations of the ground state if anhar- 
monicity is included, like in the case studied in section II. We illustrate this property below 
by studying one particular example. 
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Using Eq. ()17p . the energies of the different levels can be determined at first order in 
perturbation theory, omitting the two uniform shifts: in^^uy = ^n^,ny + ^nS^)- We consider 
a system with four electrons, two spin up and two spin down, so that only the two occupied 
and the first unoccupied levels should be considered. We first introduce a parameter, 5, 
which control the asymmetry of the harmonic confinement: uo^ = uj — 5 and uOy = uj + 5. 
The energies of the three lowest levels including the anharmonic corrections (fTTj) are found 
to be 



hfl{E) = h{2uj - 5) + f ^ + 9A§ (18) 
eo,i(^) = /i(2^ + 5) + |^ + 3A§ 

where we have introduced xe = — the shift of the center of mass induced by the electric 
field.— At a critical field, E'c, there is a change of the ground state electronic configuration, 
nny%"^' from 0^ to 0^1^. The critical field is such that l\i^{E^ = eo,i(-Ec); within our 
approximation it is given by 



^ ^ hu^ 2h5f3t - 3A 

|g| V 6A ■ ^ ^ 

This change abruptly modifies the electron density in both direction, x and y, which can be 
calculated, at our level of approximation, by using the wave functions (fT^ and (fT^. 

As already discussed for the case of hard- wall boundary confinement, for E > Ec the 
electronic cloud suddenly swells in the y-direction, perpendicularly to the applied field. The 
change of density in the transverse direction is given by 



'2 



Sp'^iy, E) = U ^e-^ly\2(3y - l)Q{E - E,) (20) 



2 V TT ' 

where Q{x) is the usual step function 

, ^ f 0, i/ X < 
e(x)= . (21) 

[l, if x>0 

Results for our example are shown in Fig. 5: electrons are transfered from the center to the 
edges of the rectangular dot, parallel to the direction of the field. 

The center of mass of the electronic system is continuously shifted by in the longitu- 
dinal direction (the x-direction) when the electric field is tuned (Fig. 4). In addition, at the 
transition, E = E^, the internal structure of the electronic cloud is modified in this direction 



too. To characterize the changes, we may consider the difference between the longitudinal 
electron density i.e. the electron density where the y-coordinate has been integrated out, 
taken at two different values of the field, and E~ , slightly above and below the transition, 
Ec, respectively: 



where p^(x) = dyp{x,y), p the electron density. Results for our example are shown 
in Fig. 6 where we have taken (3xXe = 0.5: electrons are transfered from the edges to the 
center of the dot. In conclusion, at the transition, the electronic cloud grows in the transverse 
direction (y— direction) but is slightly contracted in the longitudinal direction (x— direction). 
For Quantum Dots in GaAs, huj^ is about a few meV that corresponds to an effective Bohr 
radius, l/Px, of a few tens of nanometers; with these conditions a typical value for the 
critical field, at which the transition should occur, can be estimated with Eq. (|T!I|l . to be 
Ec ~ 10~^V.m~^ , a value commonly used in most experimental situations. 

Before closing the section, one should say that in order to observe this jump of electron 
density, the number of electrons in the dot has to be adjusted in a way to reduce the 
gap between the highest occupied level and the lowest unoccupied level; for our particular 
example, the case with 4 particles is optimal in this respect. 

IV. TWO LATERALLY COUPLED QUANTUM DOTS 

The effect pointed out above should have consequences in various experimental situations. 
In this last section we propose an experimental set-up that should enable to detect it. It is 
inspired by a proposal made in the late 90s to use semiconductor quantum dots as building 
blocks for hypothetical future quantum computers.— Following Ref. we consider two 
lateral quantum dots weakly coupled (see Fig. 6), described by the following model 



5p^(x, Et, E:) = p^(x, Et) - p^(x, E;) 




2m 



(23) 



a 



with 



Vixa,ya) = ^mujlxl + -mu;^ 




(24) 
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where the summation is over the particle indices, a and /3. The first term in brackets 
in Eq. (j^Hj) describes free electrons in a confined potential; in particular the term V{x^ y) 
represents a typical double well, which is well approximated by an harmonic potential when 
the distance 2a between the dots is large compared to the effective Bohr radius (iy Then, 



with respect to Ref. Il3| we add here the anharmonic contribution and an applied electric 
field perpendicular to the line joining the centers of the two quantum dots. The second 
term of Eq. is the Coulomb interaction supposed to be long-ranged; k is the dielectric 
constant and Tq, is the vector coordinates of the electron a, in the 3D space. 

We assume that the two dots are weakly coupled, meaning that the distance between 
them, 2a, is large (2a > 1//5); in this case, the tunneling probability of an electron from 
one dot to the other is small and the two dots can be considered, at first approximation, as 
being uncoupled. The analysis that we have done in the previous section for one single dot 
is then relevant for the coupled system. We consider the case with a small odd number of 
particles in each dot, leaving one unpaired electron in dot 1 and dot 2. In this case, under 
our assumptions, the low-energy properties of this system can be described by an effective 
Heisenberg Hamiltonian for spin 1/2 ,^1-^^1^^ 

Hueis = </Si.S2 (25) 

where the spin operators Si and S2 are associated to the spin of the unpaired electrons 
localized in dot 1 and dot 2, respectively. The effective exchange integral, J, can be estimated 
starting from the more microscopic model ()23|): it is defined as the energy difference between 
the lowest triplet and singlet states. Its magnitude depends on the Coulomb potential and 
the interdot tunneling probability. Our basic idea may be summarized as follows. We have 
seen, in the previous section, that an applied electric field along the x-direction may change 
abruptly the electron density in the y-direction: in each dot, at a critical field, Ec, electrons 
are transfered from the center to the boundaries. This phenomenon should contribute to 
increase the overlap between the two electronic clouds centered in dot 1 and dot 2, and 
therefore it should modify the magnitude of the effective exchange integral of Eq. ()25|). 
Consequently, depending on the magnitude of the electric field one expects the system to 
be described by two different Heisenberg models, with different exchange integrals, J^, for 
E < Ec, and J^, for E > E^. In other words, the exchange interaction of the effective spin 
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model is a function of the electric field that can be written as (see Fig. 9) 

J{E) = J<e{E,- E) + J>e{E - E,). (26) 

In the following, we determine these exchange integrals using a generalization of the Heitler- 
London method for the hydrogen molecule, to a number of electrons larger than twc-^; this 
type of calculation is valid for large interdot distances. Within this approximation we first 
consider the solutions for the isolated non-interacting dots (section III) before using them 
to build the proper combinations of Slater determinants to describe the whole system with 
interaction. 

The ground state of each isolated dot has been studied in the previous section, at first or- 
der in perturbation of A. We start by introducing new fermionic operators, L]^^^^^^^ {Rli^,ny,a) 
which create one electron with spin a (| or J.) in the orbital state labeled by the quan- 
tum numbers Ux and ny in the left (right) dot; they are associated with the wavefunctions 
'fn^{Px{x-XE))^ny{Py{y + a')) {^PnAPx{x - X E))^ny{Py{y - 0-)) and energies e„^,„^ . We specify 
our study to the case with three electrons in each dot with the same number of spin up and 
down {Sz = 0). In this case, it is sufficient to consider only the three lowest levels (fTHj) . and 
the gap between occupied and unoccupied states is reduced. Using the previous results, one 
concludes that the electronic configuration of the ground state changes from 0^, for E < Ec, 
to 0^1^, for E > Ec- This is depicted schematically in Fig 8, where only one particular 
configuration is shown. Starting from these solutions for isolated dots, we build for each of 
these two cases, the singlet and triplet states of the whole system. For E < Ec, we get 

n = ^{ni±n) (27) 

with 



(28) 



~ -^o,o,i-^o,o,T-^o,o,i-^o,o,T-^i,o,i-^i,o,Tl^ ^ 
^ -^o,o,i-^o,o,T-^o,o,j,-^o,o,T-^i,o,j,-^i,o,T|0 > 
and in the same way, we get for E > E^ 

with 



iTr> rt rt pt pt pt rt In ^ 

^Ti ~ -^o,o,i o,o,T-"'0,o,i-"'0,o,T-"'0,i,i 0,1,?!^ ^ 

\T/> — t"^ R'l' R'l' in^ 

^J-T ~ -^o,oa-^o,o,T-^o,oa-^o,o,T-^o,ia-^o,i,Tl"^ ^ 

12 



(30) 



where the + sign is for the singlet, the — being for the triplet, and |0 > is the vacuum. 
Next, after evaluating the triplet and singlet energy for the two different cases. 



.<,> 



<Vl/<'> 




<M;<'> 


M/<.>> 



(31) 



we can determine the two different exchange integrals, and of the effective Heisenberg 
model ^ 



J< = e< - e- 
J> = e> - 

j;iven by the difference between the triplet and the singlet energy.— 

The wavefunctions used to describe the systems are not all orthogonal since 



(32) 



(33) 



This property makes the calculation of the energy (j3ip more cumbersome. However, using 
the wavefunction (jl3|) . we can estimate that the overlap behaves as Sny,n' = (Polynomial 



of a) X e'^y""^ . Since we assume a to be relatively large with respect to 1/Py, a reasonably 
good approximation consists to keep in the calculation, terms of the type 5"™ „/ with m < 2. 
With this simplification, after lengthy but straightforward calculations, we obtain 



J< = e-^X[^i?<K( ^ 



^) + :^i?<E(- 



with 



Re 



2n 12/3^+/3^- 15/3g/3g 



and 



J> = e-'^>' [^i?>K(^^5IE) + ^ /?>E( 



^ 2(/32-/32)2 



)+-;zu;,(l-3/3y)] (34) 



(35) 



)+^.;,(2/3ja2-l)(-^ + ^/3>^-9/3^^ 



with 



R\ 



R% 



2n -/3^ + (48/34a4-48/3ga2+15)/34 + (8/3ga2+2)/3g/3g 



2/32(/32-/32) 



a 
(36) 



(37) 



2n /3g + (24/3ga4_20/3ga2+5)/3g 
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where K and E are the complete elhptic integrals of the first and second kind, respectively, 
defined as^i 

K(k) = / , E(k) = / deVl-k^ sin^e. (38) 

Jo ^/l-k^sm^e Jo 

Although, the two expressions ()34|1 and ()3fi|l are only valid in the regime of weak anhar- 
monicity, A//?^ << hu^/y, and large interdot distances, Pya > 1, they allow us to illustrate 
our purpose. Indeed, they show that the exchange integral can be significantly modified by 
applying an electric field perpendicular to the axis joining the two dots. Taking relevant 
numerical values for dots in GaAs,^^ huy = 3meV, which corresponds to l/(3y = 20nm, 
and (3yq^ / n/Tiujy = 2.1 and specifying our system with Pya = 1.3, we get for Px/Py = 0.9, 
J< = O.lSmeV and J> = 0.6meV and, for P^/Py = 0.95, J< = 0.19meV and J> = 3.13meV. 
We estimate the critical field in the latter case by using Eq. (fnijl : assuming X/P^ = O.lSmeV, 
we obtain = 35.10^''Vm^^. In both cases, the energy difference between the singlet and 
triplet states is large enough to be probed experimentally, as it was done in Ref. |l5| for two 
coupled dots with an exchange energy tuned by an applied magnetic field as suggested in 
Ref.ll^. 

D. Loss and D. P. Di Vincenzo proposed to use a system of two weakly coupled quantum 
dots as possible realization of quantum gate.-^ With one unpaired electron in each dot, the 
system is described by an effective Heisenberg model such as (j^ . By changing the value 
of the exchange energy with external sources in such a way that JJ^ dtJ{t)/h = vr, it would 
be possible to perform a swap operation that exchanges the quantum state of dot 1 and 
dot 2. When combined with single-qubit operations, this operation can be used to build a 
quantum XOR gate that was earlier recognized to be a universal quantum gate^^: a XOR 
along with single-qubit operations may be assembled to do any quantum computation. Thus, 
the study of quantum gates is reduced to the study of the exchange integral of the two dots 
system and in the way it can be controlled experimentally. Two main approaches have been 
proposed up to now to control J{t): (i) by a local magnetic field perpendicular to the plane 
of the dots and, (ii) by a local electric field parallel to the axis joining the center of the 
two dots. Our proposal offers another possibility. With respect to the method described 
in Ref. 3, the main advantages of the mechanism proposed here are (i) to let the average 
charges unchanged within the two dots, quantity that would be modified by an electric field 
applied along the y-axis (see Fig. 7); (ii) it is easier to apply a local electric field than a local 
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magnetic field. 

As a remark, we would like to point out that there is no need to consider identical dots. 
Both the number of particles within each dot and the geometry of the dots might be different. 
In such situation, the exchange integral, J{E), should show two discontinuities instead of 
one at two different critical field, E^i and Ec2, characterizing the electronic changes in the 
dot 1 and dot 2, respectively. 



V. CONCLUSION 



To conclude, we summarize our main results. We have considered systems of free-electrons 
trapped in a small region of space. Two confining potentials were studied as examples: hard 
wall boundary conditions and harmonic potential with an additional perturbative term. 
In both cases, an electric field was shown to cause abrupt changes in the electron density at 
a critical value, Ec- These changes occur not only in the direction parallel to the field but 
also in the direction perpendicular to it; they can be related to finite size effects. Exploiting 
this outcome, we have then considered a system of two weakly coupled Quantum Dots and 
shown that an applied homogeneous electric field is able to change abruptly the difference 
in energy between the lowest singlet and triplet states. This property may be important 
in the context of quantum computation, since it could be exploited to realize a quantum 
gate. Our results were obtained by using deliberately oversimplified models, our aim being 
to point out basic mechanisms that should survive more realistic descriptions. Of course, to 
be more quantitative, more detailed models should be studied in the future. 

As a last remark, we have seen that the electric field induces degeneracy in the electronic 
spectrum at E = Ec- In this case, the Coulomb interaction should play a central role that 
could significantly modify our conclusions for fields close to the critical value: in particular, 
the total spin of the system may be changed. This kind of transitions was already studied, 
for instance, in Ref. [isl, for cylinders pierced by an Aharonov-Bohm flux, and in Ref. for 
single-wall carbon nanotubes with an inhomogeneous electric field. 
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FIG. 1: Rectangular potential well with infinite walls [L^ > Ly). A uniform electric field, E, is 
applied along the long axis. 
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FIG. 2: The two lowest bands, ny = 1 and Uy = 2, of a system of free-electrons confined in a 
rectangular potential well such that Lx = lOLy. The long dashed curves are for E = 0, the full 
curves ior E = The positions of the bands with electric field are shifted down by tt^, for 

clarity. In both cases, the dots denote the position of the discrete set of quantum states available 
due to finite sizes. They are labeled by the integer n, equivalent to Ux in the case without electric 
field. The open diamonds and the open circles are the doubly occupied levels in the case with 
34 electrons and Sz = 0, without and with electric field, respectively. The electric field induces a 
change in the set of occupation numbers, the ground state changing from 1^^ to 12^2^ (see text): 
two electrons jump from the band nj, = 1 to the band Uy = 2. 
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FIG. 3: Changes in the transverse electron density, 6p^{y,E > Ec), of a system of free-electrons 
confined in a rectangular box with hard-wall boundary conditions. At a critical value, Ec, an 
electric field induces abrupt changes in the set of occupation numbers, An^, defining the ground 
state (see text) and, consequently, in the electron density in the direction perpendicular to the 
applied field. In this example, the ground state changes from 1^^ to 1^^2^. The value of the critical 
field depends on the geometry of the well: it is Ec ^■^ 2m\q\ ~ ^^^y -^c — 0-01 2m\q\ 

for Lx = 9.8Ly, for instance. 
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FIG. 4: Rectangular quantum dot with harmonic confinement potential, such that uix < (^y (full 
line). An applied electric field along the x-axis shifts the position of the center of mass of the 
electronic system hy xe = (dashed line). 
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FIG. 5: Changes in the transverse electron density, 5p^{y,E > Ec), of free-electrons trapped in a 
harmonic potential with an additional weak anharmonic term, Xx'^. In our example, the ground 
state configuration, H^y"^' changes from 0^ to 0^1^, at £^ = Ec- As a consequence of that, some 
electrons leave the center of the dot towards the edges parallel to the direction of the applied field. 
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FIG. 6: Changes of electron density at E = Ec, in the longitudinal direction (parallel to the applied 
field) of free-electrons trapped in a harmonic potential with an additional weak anharmonic term, 
Ax^. In our example, the ground state configuration, 11%"^' changes from 0^ to 0^1^, at = Ec- 
The dashed lines show the longitudinal density at a value of the field slightly smaller than Ec- 
The dotted-dashed lines show the longitudinal density at a value of the field slightly higher than 
Ec- The full lines show the changes in the longitudinal density at the transition: electrons are 
transfered from the edges to the center of the dot. Here, we take PxXe = 0.5 (see text). 
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FIG. 7: Two laterally coupled rectangular quantum dots with a harmonic confinement such that 
ujx < (full lines). An applied electric field along the x-axis shifts the position of the center of 
mass of the two quantum dots, supposed to be identical, hy xe = (dashed lines). 
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FIG. 8: Example of electronic configuration of two weakly-coupled quantum dots with three elec- 
trons in each. An applied electric field can change the energy order of the one-electron states, 
{nx,ny), in each dot. a) For E < Ec, the states (0,0) and (1,0) are populated, b) For E > Ec, the 
states (0, 0) and (0, 1) are populated. At the critical field, Ec, the unpaired electrons are transfered 
from the state (1, 0) to the state (0, 1) in each dot. 
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FIG. 9: Exchange integral of the effective Heisenberg model (|25|) . describing the low-electronic 
properties of two weakly coupled quantum dots with one unpaired electron in each, as function of 
the applied electric field. At a critical field, Ec, the electronic structure of each dot is modified 
which causes a sudden jump of J, changing from to J^. This property may be used to perform 
a swap operation (see text) in the context of quantum computation. 
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